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ABSTRACT 

Context. The solar surface magnetic field is connected with and even controls most of the solar activity phenomena. Zeeman effect 
diagnostics allow for measuring only a small fraction of the fractal-like structured magnetic field. The remaining hidden magnetic 
fields can only be accessed with the Hanle effect. 

Aims. Molecular lines are very convenient for applying the Hanle effect diagnostics thanks to the broad range of magnetic sensitivities 
in a narrow spectral region. With the UV version of the Zurich Imaging Polarimeter ZIMPOL II installed at the 45 cm telescope of 
the Istituto Ricerche Solari Locarno (IRSOL), we simultaneously observed intensity and linear polarization center-to-limb variations 
in two spectral regions containing the (0,0) and (1,1) bandheads of the CN E'l, - X^l, system. Here we present an analysis of these 
observations. 

Mettiods. We have implemented coherent scattering in molecular lines into a NLTE radiative transfer code. A two-step approach was 
used. First, we separately solved the statistical equilibrium equations and compute opacities and intensity while neglecting polariza- 
tion. Then we used these quantities as input for calculating scattering polarization and the Hanle effect. 

Results. We have found that it is impossible to fit the intensity and polarization simultaneously at different limb angles in the frame- 
work of standard ID modeling. The atmosphere models that provide correct intensity center-to-limb variations fail to fit linear polar- 
ization center-to-limb variations due to lacking radiation field anisotropy. We had to increase the anisotropy by means of a specially 
introduced free parameter This allows us to successfully interpret our observations. We discuss possible reasons for underestimating 
the anisotropy in the ID modeling. 
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1. Introduction 

The Q/I spectrum formed by coherent scattering is a rich source 
of information about the solar atmosphere and was called there- 
fore the "second solar spectrum" (Ivanov 1991 1. The proper- 



ties of scattering processes are modified in a magnetic field via 
the Hanle effect (Stenfio 1982), and it makes the second solar 
spectrum a unique tool to access weak, spatially unresolved en- 
tangled magnetic fields with mixed polaiity, which cover about 
99% of the photospheric volume and cannot be measured with 



Zeeman effect because of signal cancellation (Stenflo 2004 
[Trujillo Bueno et al. 2004| ). 

Scattering in molecular lines plays an important role in form- 
ing the second solar spectrum and even dominates in some spec- 
tral regions. The sensitivity of molecular lines to the Hanle effect 
varies significantly with the total angular momentum of the up- 
per level J. Since a narrow spectral region can contain lines with 
quite different J values, this allows us to employ the differen 
tial Hanle effect (cf |Stenflo et al.|[T998t |Trujillo Bueno||2003 
[Berdyugina & Ruri|2004| i 



Shapiro et al. ( 2007b[ l developed the theory of the Hanle ef- 
fect in the CN violet (B^l, - X ^2) system, which is a signif- 
icant feature of the solar spectrum, including Q/I, and is very 



convenient for applying the differential Hanle effect diagnos- 
tics. Then, while employing a simple radiative transfer model, 
Shapiro et al. ( 2007a| l applied this theory for modeling CN vio- 
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let system lines in the second solar spectrum and for determining 
the turbulent magnetic field strength. Despite the successful fit in 
several spectral regions, this simple model failed to fit strong CN 
lines. Moreover, it was impossible to self-consistently interpret 
the differences in the deduced model parameters obtained from 
different spectral regions. To solve these problems, we present 
here a new, more self-consistent model of ID polarized radiative 
transfer. 

For our purpose, we carried out simultaneous observations 
of center-to-limb variations in the intensity and linear polariza- 
tion in the CN (0,0) and (1,1) bands. Such observations are very 
useful for constraining the atmosphere model and studying the 
height dependence of the deduced parameters as observations at 
different Umb angles sample different heights in the solar atmo- 
sphere ( [Faurobert & A rnaud 2003). We developed a new model 
to interpret these observations. Here we present the first results 
of this interpretation and discuss the general problems of simul- 
taneous modeling of I/I^ and Q/I in ID . 

In Sect.|2]we briefly describe our new observations. Then, in 
Sect. [3] we introduce the radiative transfer model and report our 
main assumptions. In Sect. |4] we discuss the problem of lacking 
radiation field anisotropy arising in simultaneous ID modeling 
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of I lie and Qjl and show that it is not possible to fit all observed 
data in the framework of a ID one- or multi-component atmo- 
sphere model. In Sect. |5] we discuss possible physical mecha- 
nisms leading to the insufficient anisotropy in ID modeling and 
introduce a mean to increase the anisotropy and solve the prob- 
lem. Finally in Sect. |6]we show that modeling with an increased 
anisotropy can provide us with a good fit quality and allows us to 
measure the magnetic field strength. We also discuss the model 
dependence of the deduced magnetic field value. 



2. Observations 

With the UV version of the Zurich Imaging Polarimeter 
ZIMPOL II installed at the 45 cm telescope of the Istituto 
Ricerche Solari Locarno (IRSOL) we simultaneously observed 
Stokes / and Qjl center-to-limb variations in two wavelength re- 
gions, which contain the (0,0) and (1,1) bandheads. Observations 
of these two regions were made on December 13 and December 
14, 2007, respectively. For all observations the slit subtended 3.4 
arcmin and was placed parallel to the nearest limb. The slit width 
corresponded to a spatial resolution of 1 arcsec. Exposure time 
for one exposure frame was 5 sec, but to reach the required num- 
ber of photons in order to get the desired signal to noise, the total 
exposure time at different /i positions was varying from 20 min- 
utes to 48 minutes. Such a choice of the wavelength regions is 
justified for magnetic field diagnostics as bandhead regions con- 
tain a mixture of lines with different /-numbers and accordingly 
different magnetic sensitivity and formation height. Therefore, 
the differential Hanle effect technique can be applied. Moreover 
the bandhead is formed higher in the atmosphere than the radia- 
tion at the nearest wavelengths (1-2 A away from the bandhead) 
so the information about conditions at a broad range of heights 
can be inferred. 

The intensity and linear polarization in the region containing 
the (1,1) bandhead (3868.9-3871.4 A) was observed at five limb 
angle positions from ^ - Q.\ lo ^ - 0.5 (here ^ is the cosine of 
the angle between the propagation direction of radiation and the 
local solar radius). The intensity in the second region containing 
the (0,0) bandhead (3879.7-3883.6 A) was observed at ten Hmb 
angle positions from // = 0.1 to // = 1.0 with the constant step 
Aji The polarization in the second region was observed at 

two positions: // = 0.1 and ji - 0.2. The pointing of the telescope 
allows a precision of about 1 arcsec, but due to seeing and the 
curvature of the solar limb relative to the straight spectrograph 
slit, we need to take into account an uncertainty of about ±0.02 
si H = 0.1. The error for jj. values quickly decreases toward the 
disk centre. 

The zero level of the polarization scale is uncertain in the ob- 
servations. In Sect. 3.2.2|we present an algorithm which allows 



us to simultaneously determine the continuum polarization and 
the zero level of polarization scale. 



3. Radiative transfer model 

In the following we introduce our numerical method employed 
for solving the radiative transfer problem. We have performed 
calculations in a framework of plane-parallel model atmo- 
spheres. The modeling consists from two main steps, as was 
previously suggested by |Fluri et al. ( |2003l ). We firstly calculate 
opacities and intensity without taking into account any polar- 
ization. Then we iteratively calculate polarization assuming that 
opacities obtained in the first step remain unchanged. 



3.1. Phase matrix 

In the weak-field regime the redistribution matrix can be factor- 
ized into a scalar function depending only on frequencies of in- 
coming and outgoing radiation and the phase matrix. The phase 
matrix does not directly depend on frequencies, but it has differ 
ent appearance in different frequency domains ( jBomimer|1997 
IFluri et al.|20"03| ). 

Under the assumption of complete frequency redistribution 
(CRD) and an isotropic single-value turbulent magnetic field the 
phase matrix for the line core domain (roughly speaking in this 
domain both absorption and emission of a photon occur close to 
the line core) is given by ( Fluri et al.|2003| l 



P(M,I^', W2, Wh) = En+ 3/4 W2WH ■ P' 
and in the line wing domain by 

P(M,I^', W2) = £11 + 3/4 W2 ■ P'^^'K 



.(2) 



(1) 



(2) 



In these equations fi and /i' indicate the outgoing and incoming 
directions of the scattered photon, W2 is the effective scattering 
polarizability, Wh is the Hanle depolarization factor, and 
and f'^* are 2 X 2 matrices (as we solve the radiative transfer 
problem for Stokes / and Q only) which are given by 



En = 



and 



1 




(3) 



P''\^i,^^') = 



1 / kl - 3//2)(l - 3/2) (l-3//2)(l-//'2)' 



2\ (l-/)(l-3/0 3(l-/)(l-//'0 



(4) 



In agreement with the principle of spectroscopic stability, 
line wings are not affected by the Hanle effect, and the Hanle 
depolarization factor Wh does not enter Eq. Values of the 
dimensionless frequency for a core- wing boundary separation 
are calculated for certain values of the Voigt parameter a by 
jBommier ( 1997[ l. For a given damping parameter a, we inter- 
polate between the known values. 

The matrix Eu represents isotropic, unpolarized scattering, 
while /"'^^ describes linearly polarized coherent scattering (e.g., 
Stenflo|1994| l. Due to the latter even initially unpolarized radia- 
tion becomes polarized after scattering, and the degree of polar- 
ization scales with the effective polarizability W2. 

The Hanle effect in the line core for the case of a turbulent 
single-value magnetic field is introduced by ( Stenflo|I982 1 



Wh = 1 - 0.4 



where 



Th 



0.88 





l+ri"l+4y2 



5lS 



Fr + Fi + flf ■ Fe 



(5) 



(6) 



Here Fr,Fi and Fe are the radiative damping, inelastic and elas- 
tic collision rates, respectively, d is the efficiency of the depolar- 
ization by elastic collisions, B is the magnetic field strength, and 
is the effective Lande factor of the transition's upper state. 
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3.2. Source function 

The total source function S = S(t, v, fj) consists from parts aris- 
ing due to thermal and scattering processes in the continuum and 
spectral lines. It is given by 



S = 



n 

1=1 



!=1 



(7) 



where ffc = fcc + cTc is the continuum opacity (here k^^ and ct;, 
are the continuum absorption and scattering coefficients), a'^ - 
kg + cr'f is the total line opacity in the i-th line (here k'^ and cr'^ are 
the line absorption and scattering coefficients), and ip'(v) is the 
normalized Voigt profile function. Sc^dT, v,fi) is the continuum 
source function, and 5^ = 5'^(t, v,//) is the line source function of 
the i-th line. 



collisions or by the absorption of the photon in any j-th line line 
which shares the upper level with the i-th line (see Fig.[Ti. Each 
excitation possibility gives rise to a separate term in the line 
source function. So the source function in the i-th line within 
the CRD approximation can be written as 



S'( - ^ Ssc + ^Ih-^th' 



(10) 



where s'sl and are the branching coefficients (which are depth 
dependent) and 5^'^ is the scattering integral which corresponds 
to the absorption in the j-th line and emission in the i-th line. It 
is given by 



+ 00 +1 



P(ii,lj', W'^\ W'ff)I(v',iu')(p'(v') d/u'dv'. 



(11) 



3.2.1. Continuum source function 

The continuum source function is given by 



In Eq. (10 1 scattering integrals 5^^ with / 9^ j describe the 



(8) 



Here the vector Bth - (fiv, 0) is the thermal source term, which 
is given by the Planck function. The monochromatic scattering 
in the continuum is introduced by 



1 



(9) 



where / is the Stokes vector of the radiation field, and P is the 
phase matrix, which is given by Eq. (|2j. The coefficient Wf^ is 
the effective continuum polarizability. 

According to Eqs. (j2]i and (j9]l the continuum polarization is 
proportional to the effective continuum polarizability. The exact 
value of the latter depends on many factors which are model de- 
pendent (see |Stenflo|2005 1. The situation is further complicated 
by the line blanketing effect. To modify the continuum absorp- 
tion coefficient k^ we used the value of the opacity correction 
factor from Busa et al. ( 2001] l, which was adjusted to better fit 
the intensities. Due to the lack of additional information we em- 
ployed the same value to modify the continuum scattering coeffi- 
cient CT;,, which is not necessarily correct. The possible error will 
alter the scale of the continuum polarization. Therefore, the con- 
tinuum effective scattering polarizability was considered as 
a free parameter noting however that it also contains a possible 
error from the line blanketing coefficient. 

The value is determined by fitting the observed center- 
to-limb variations of the continuum polarization. However to ob- 
tain the absolute value of polarization, the observed Qjl profiles 
have to be shifted by the value of the zero level of polarization 
scale (see Sect.[2|. We choose the shift for the observed QII and 
the W^** value so that the observed continuum polarization level 
coincides with the calculated. While the value of the shift is the 
same for all //-values, the influence of the effective polarizability 
on the continuum polarization depends on /i. As we have obser- 
vations at several yu-values, these parameters can be determined. 
The obtained value of the effective continuum polarizabiUty is 
about Wf 



0.6 for both spectral regions. 



3.2.2. Line source function 

For the photon to be emitted in the i-th line, the upper level of 
this line has to be excited. The excitation can be caused by the 



Raman scattering, while the scattering integral S^^. describes the 
Rayleigh scattering. The last term corresponds to the thermal 
photons, emitted after the collisional excitation. 

To illustrate the connection between the anisotropy of the 
radiation field and scattering polarization it is useful to express 
the scattering integrals in the form 



5- = Jr'^(y' 



)ip'(v')dv'. 





where 



+ 1 



(12) 



P{M,l.^',w'2',W'H)I(v',i^')dfi'. (13) 



Substituting the phase matrix from the Eq. ([Tl-Q we obtain 
V2 ■ ■ 

e;\v') = y»(v') + — w^'-'Wh(3//2 - 1) ■ JliV), 



where J^iv") is the mean intensity 

+1 



(14) 
(15) 

(16) 



and J^iv') is the second moment 

+1 

4(^')=j^ J((3yu" - l)/(v',//')+30u'' - l)Q{v',fi'))dfi'. (17) 

-1 

The first term in Eq. ( [T4| describes isotropic scattering and 
depends only on the mean intensity. The second term repre- 
sents coherent scattering and depends on the magnetic field. It's 
value is negligible in comparison with the first term, because the 
anisotropy of the radiation field is usually very small. Therefore, 
Stokes / is mainly determined by the mean intensity. In contrast, 
the polarization part of the line source function (see Eq. (TSj)) 
does not contain contributions from isotropic scattering and lin- 
early depends on J^y Due to the Wh factor this part is also subject 
to the Hanle effect. 
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Therefore, the (2//-signal is proportional to the anisotropy of 
the radiation field which is given by (see Trujillo Bueno|200r 
IHolzreuteret al.|2005| l 



A - 



J 



(18) 







By definition, the anisotropy describes the excess of the radiation 
incident from the radial directions (close to vertical) in compar- 
ison with the radiation incident from the side (horizontal direc- 
tion). 



3.3. Two-step approach 

For solving the statistical equilibrium and the radiative trans- 
fer equations for the CN violet system lines we employ the nu- 
merical NLTE code (hereafter RH-code) written by )Uitenbroek] 
( |2001| l and based on the multilevel accelerated lambda itera- 
tions (MALI) method (cf. |Rybicki & Hummer] [T99 1 , ,1992| ). 
Contributions from other lines were considered as "background" 
and treated in LTE. The statistical equilibrium equations in 
the RH-code are solved under the assumption of LTE pop- 
ulations within a single vibrational level since it would be 
time-consuming to solve them for each individual vibrational- 
rotational level (as we have to deal with more than 2000 transi- 
tions simultaneously). Such an approximation is commonly used 
in NLTE molecular calculations (cf. [Thompson 19731 IMount 



& Linsky|1975[|Ayres & Wiedemann! 1 989[ |Uitenbroek[2000| l, 
because oscillator strengths of pure rotational radiative transi- 
tions are negligibly small and radiative processes can contribute 
to the rotational population balance only via two-step processes 
(Raman scattering for the ground electronic state and emission 
followed by absorption for the exited state). 

We have included calculations of electronic-vibrational 
molecular transitions with the fine structure into the RH-code in 
order to use it for the CN violet system computation. Although 
the modified code can be used for a more general case, for the 
CN violet system we consider only diagonal vibrational bands as 
the Franck-Condon factors of the non-diagonal bands are much 
lower than those of the diagonal. 

With the RH-code we compute the opacities and intensity, 
neglecting polarization. Then, we use these as the input for the 
second code (hereafter POLY-code) written by [Fluri & Stenflo 



( |200 3 ) and Fluri et al. (2003), which iteratively solves the polar- 
ized radiative transfer equation, taking into account the Hanle ef- 
fect and assuming that opacities obtained in the RH-code remain 
unchanged (which is a good approximation since the degree of 
polarization in our calculations is always lower than 1%). The 
POLY-code was adjusted for dealing with many blended lines as 
initially it was designed for treating only a few non-overlaping 
atomic lines, while in molecular bands even a narrow spectral 
region can contain several hundred lines, which have to be com- 
puted simultaneously. 



3.4. Thermal-scattering branching of the source function 

As in most NLTE-codes, the line source function in the RH-code 
is calculated as a function of NLTE populations, temperature 
and frequency (see Uitenbroek| [2000 ). The source function de- 
viations from the LTE are defined via population departure coef- 
ficients (the ratio between the LTE and NLTE populations). Such 
a formalism is sufficient for calculations of the intensity field but 
does not provide any direct information about the contribution 



N+ 1 



N-1 




Fig. 1. Scheme of the radiative excitation between the fine struc- 
ture levels of the CN violet system. Doublet states are marked 
with the quantum number A^. Allowed transitions in the Pi, Ri, 
and ^2i2 branches are indicated. 



of scattering processes to the line. However, for polarization cal- 
culations it is important to know the exact balance between scat- 
tering and thermal processes, since polarization is produced only 
by scattering. Therefore in the POLY-code the source function is 
calculated as a sum of scattering and thermal parts (see Eq. 



10 1) 



The conversion from the RH to the POLY-approach is carried 
out by choosing branching coefficients e'^^ and ejj^ in Eq. JlOb. 
Defining these coefficients is a standard problem for a two- 
level system, which is described in many textbooks (e.g.'Mihalasj 
1970). However, in the case of the CN violet system the situation 
becomes much more complicated as we have to consider a lot of 
levels, many of which are collisionally and radiatively coupled 
with each other 

To find the branching coefficients let us consider the pre- 
history of the photon which was emitted in the i-th line, when 
molecule changed from the upper to the lower state. If the ex- 
cited state was collisionaly populated (it could be both the col- 
lisional excitation from a lower level or de-excitation from a 
higher energy level) a thermal photon, which contributes to the 
last term of Eq. 



10 1, will be emitted. However, if it was radia- 



tively populated by the absorption of the photon in the j-th line 
(which shares the upper level with the i-th line), a nonthermal 
photon, which contributes to the s'sc Sfj term in Eq. J 10b, will 
be emitted. Therefore, if Ctot and /?tot = Z / are the total rates 
of the collisional and radiative processes populating the upper 
level, then the part of thermal photons in the i-th line is given by 



^th = Ctot /(Ctot + ^tot), 

while the other part will be nonthermal: 

<5sc = ^tot/(Ctot + ^tot)- 



(19) 



(20) 



As we solve the statistical equilibrium equations only for vi- 
brational levels, the coefficients 6th and 6sc are also calculated 
only for vibrational levels of the excited electronic state B-E 
(adding up all significant radiative and collisional processes, 
populating the considered vibrational level). We assume that the 
same fractions of thermal and scattered photons 6±(v) and (Jsc(v) 
apply to every rotational level within a given rotational band 
(v,v). This seems to be a good approximation as the rotational 
energy is low in comparison with the electronic energy of the 
excited B^S state. 

In the case of a two-level system the branching coefficients 
ejjj and s'^l. are equal to the coefficients 6th and 6^,^. However in 
our case of a multi-level system this equality is not valid and 
the branching coefficients depend not only on the coefficients 
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Fig. 2. Observations (thick solid lines) of Stokes ///c and Qjl in the region of the (1, 1) bandhead and fits with the FALC atmosphere 
model. The assumed magnetic field strength in the calculated spectra is G. Two curves presented in each picture correspond to the 
NLTE (thin solid) and LTE (dotted) calculations. 



5th and (5sc> but rather on the whole balance of populating and 
depopulating rates. 



levels from RH, we firstly calculate the line source function with 
the standard equation 



hJ 



The following algorithm was employed to compute the e! 
and ej^ coefficients. Using the NLTE populations of the and w, 



n,,A 



Ujlj 



(21) 
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where A 1,1' and Bw are the Einstein coefficients. The source 



function given by Eq. (21 1 has to be equal to the intensity part of 
the source function given by Eq. ([TO]). So one can write 



',Bu 



(22) 



where 7?*^ is the intensity part of the scattering integral from 
Eq. ([11).' 

Secondly, the ratio between the thermal (arising from the col- 
lisional excitation) and j-th nonthermal (arising from the radia- 
tive absorption in j-th line) contributions to the source function 
is given by 



'Jrsc 



Ctot 



(23) 



The Eq. ( 23 1 can be written for each from the nonthermal con- 
tributions to the source function and therefore, together with the 
Eq. ( 22 1 the number of the equations for the branching coeffi- 



cients is equal to the number of the branching coefficients. Since 
the equations are not linearly dependent the branching coeffi- 
cients can be uniquely determined. We note that in our case of 
the radiative transfer problem with the multi-level coupling the 
sum of the branching coefficients is not necessarily one (see the 
lower panel of Fig. 3). 

For every line in the CN violet system there are three differ- 
ent ways to radiatively populate the upper state before the emis- 
sion (see Fig.[T]l. One of them, in a main branch, corresponds to 
Rayleigh scattering and the two others to Raman scattering (one 
in a main branch and one in a satellite branch). The oscillator 
strength of the transitions in a satellite branch are approximately 
four orders of magnitude weaker than the transitions in the main 
branches and therefore the satellite transitions can be neglected. 
So in our calculations the line source function consists from two 
scatterings terms (one corresponding to the Rayleigh and another 
to the Raman scattering) and one thermal term (see Eq. ([TOjl and 
the second equality in Eq. ( 22 1). 



where f5 - AE/kT is the excitation parameter. 

As there are several collisional agents (most importantly 
neutral hydrogen and electrons) the collisional rates depend on 
many free parameters. However for computing the branching co- 
efficient and scattering polarization only the total collisional rate 
(independently on collisional agent and type) is relevant (see 
Sect. 13.41). We assume then that the total collisional rate can be 



also calculated with Eqs. (24 1-(25 1, using the partial pressure of 
neutral hydrogen (as the collisions with neutral hydrogen seems 
to be the strongest and moreover the depth dependence of the 
electron particle density is close to the one of neutral hydrogen 
in the region of the CN lines formation). 

With such an approximation the collisional rate depends only 
on two free parameters: A and B, which define the tempera- 
ture and density dependence, respectively. At the CN line for- 
mation height the density drops much faster than the tempera- 
ture. Therefore, the parameter A only slightly affects the depen- 
dency of the collisional rates on height in the atmosphere, in 
particular also when considering the small exponent of the tem- 
perature in Eq. (|24]). We found that the best fit quality can be 
reached if we put the value of A parameter equal to zero. Let us 
notice, however, that even under this approximation the branch- 
ing coefficients Sth and 6th still depend on temperature due to the 
Einstein coefficients, which connect excitation and de-excitation 
collisional rates and affect the level populations. 

Finally we end up with only one free parameter B which de- 
fines the dependency of the collisional rate on height. For con- 
venience we have replaced it with the 6th coefficient at the tem- 
perature minimum layer (hereafter 5™"). Considering the strong 
limitations of the current collision theory with molecules, the 
unknown collision parameters for CN, and our need to know 
only the total collision rate such approximation becomes very 
practicle and improves the stability of the procedure to fit ob- 
servations. We can control the collision rate with a single free 
parameter and nonetheless account for the height dependence 
via the Landau-Teller expression. In fact, by fitting observations 
we can even gain empirical constrains for collisional rates with 
the CN molecule. 



3.5. Collisional rates 

There are three main types of inelastic collisions with molecules. 
The first one alters only the total angular momentum number J 
of the molecular state, the second also the vibrational state, and 
the third the electronic state. There are many theoretical approx- 
imations and parametrizations for the rates of these three colli- 
sion types (cf. [T hompson 1973; Hinkl e & Lambert|1975[|Ayres| 



& Wiedemann 1989). The accuracy of these approximations is 



very low, and the parameters for the CN molecule are unknown. 
Moreover it is not possible to distinguish the elastic and inelastic 
collisions from our observations. Therefore, we neglect elastic 
collisions. It does not influence our conclusions but can slightly 
affect the value of the deduced magnetic field strength. 

In our calculations we used the Landau-Teller formula for 



the relaxation time of the excited state (Ayres & Wiedemann 
[1989] ) 



\n(PxtcN-x)=AxT-''^ -B 



(24) 



where Px is the partial pressure of the collision partner "X", 
fcN-x is the relaxation time of the excited CN state, T is the 
temperature, and Ax and Bx are free parameters. The associated 
collisional de-excitation rate per one molecule is 



C = (fcN-x(l-exp(-yS))- 



(25) 



4. Center-to-limb variations with standard modeling 

4.1. One-component atmosphere modeling 

As a first step we performed computations in a frame of a 
plane-parallel one-component model atmosphere, namely model 
FALC by [Fontenla etal] ( [T993| ). In Fig. ^ both the LTE and 
NLTE calculations of the intensity and linear polarization for 
different limb angles in the region of the (1,1) bandhead are 
presented. The region contains about 100 CN lines and several 
atomic blends, among them two strong iron lines around 3869.5 
A. These blends in principle have to be treated in NLTE too, but 
for simplicity we treated all blends in LTE, which reduces the 
quality of the fit at wavelengths affected by blends but does not 
influence our conclusions. 

From Fig. [2] one can see that the calculated and observed 
linear polarization curves have completely different shapes. 
Surprisingly, the calculated polarization dramatically decreases 
and for high ju-values even becomes negative (perpendicular to 
the solar limb) in the bandhead at 3871.2 A, completely con- 
tradicting the observations. Here we have chosen the collisional 
coefficient for the NLTE calculations (5™" equal to 0.15, so that 
the overall level of polarization outside the bandhead fits the ob- 
servations. The overall polarization level can be increased by 
choosing a smaller collisional coefficient, but this will make the 
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situation at the bandhead polarization even worse as its absolute 
value will increase but it will remain negative. Even excluding 
the particular bandhead region, the calculated shape of the Qjl 
curve does not correspond to the observed one and it can not be 
corrected with the Hanle effect (see Sect.|6]l. 

While the calculated polarization shows such a strange be- 
havior, the quality of the intensity fit is good for all five yU-values. 
In the LTE calculations lines are too deep, while taking into ac- 
count the NLTE effects makes lines a bit weaker and leads to 
the good agreement with observations. It is interesting that the 
higher /i-value, the weaker NLTE corrections. This is connected 
with the fact that observations at greater limb angles sample 
deeper layers of the solar atmosphere where the density and col- 
lision rates are higher and, therefore, the deviations from LTE 
become smaller 

Thus we find that the standard FALC model is able to repro- 
duce the center-to-limb variations for the intensity but dramati- 
cally fails for polarization. The good fit to Stokes / was expected 
because the FALC model has been empirically optimized to re- 
produce the solar intensity spectrum. 

In the following paragraphs we discuss the reasons why cal- 
culations with the FALC model atmosphere yield good results 
for the intensity spectrum but not for the second solar spectrum. 
The main point to recall is the fact that intensity and scatter- 
ing polarization are influenced by different physical quantities. 
As discussed in Sect. 3.2 scattering polarization scales with the 
anisotropy of the radiation field, which has negligible impact on 
Stokes / in the case of the solar atmosphere. The discrepancy 
in the fits to the intensity and polarization spectra can thus be 
understood in the following terms: the FALC model correctly 
reproduces the mean intensity and thus Stokes / within the solar 
atmosphere, but yields improper radiation field anisotropy and 
thus very bad fits to scattering polarization. 

For a deeper understanding we briefly review the main fac- 
tors that influence the anisotropy of the radiation field within a 
stellar atmosphere. The anisotropy mainly results from two com- 
peting contributions, namely the source function gradient and 
a surface effect ( Trujillo Bueno|[2"00I[ |Holzreuter et al.|[2005| l. 
The source function gradient is, under LTE conditions, due to 
the temperature gradient within the atmosphere. The tempera- 
ture drop in the photosphere leads to limb-darkening. This means 
that at any point within the photosphere the upward moving ra- 
diation field is strongest in the vertical direction and reduces to 
more inclined directions from the side, which by definition con- 
tribute to positive anisotropy. The surface effect, on the other 
hand, is active only within one mean free photon path from the 
top of the atmosphere and affects the downward moving part of 
the radiation field. Within these topmost layers the optical path 
length between the top and some point within the atmosphere be- 
comes significantly greater for inclined directions. As a result the 
downward directed radiation field is limb-brightened and thus 
contributes to negative anisotropy simply due to the presence of 
the surface. 

The two effects lead to an anisotropy that in general depends 
on the height within the atmosphere and on wavelength because 
the source function gradient is height and wavelength dependent. 
Furthermore, the influence of the surface effect is zero deep in 
the atmosphere, reaches a maximum around an optical depth t - 
1 and reduces towards the top because the absolute amount of 
downward radiation becomes negligible very close to the top. 

Figure |3] illustrates both competing factors, their influence 
on anisotropy, and their subsequent effect on scattering polar- 
ization for the case of spectral modeling of the (1,1) band in 
the CN violet system based on the FALC atmosphere. From the 
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Fig. 3. Upper and middle panels: Normalized contribution func- 
tions of Stokes / as well as the anisotropy at 3871.2 A in the (1,1) 
bandhead (thin soHd lines), at 3869.1 A in a CN blend (dotted 
lines) and at the local continuum (dashed lines). The contribu- 
tion functions are given for five values: from 0.1 (rightmost) 
to 0.5 (leftmost). The FALC temperature as a function of height 
is given in the upper panel (thick solid line, right scale). Lower 
panel: Branching coefficients for the source function of Pi (14) 
(3869.1 A) line: thermal (solid line), Raman (dotted line) and 
Rayleigh (dashed line). Assumed 5"? 



5™" value is 0.3. 



contribution functions of Stokes / we find that our observations 
sample a wide range of the solar atmosphere from the very bot- 
tom of the photosphere up to about 600 km. This is mainly due 
to the increased opacity when approaching the bandhead, which 
contains a large number of blended CN lines and thus forms at 
higher atmospheric layers. We consider anisotropy and contribu- 
tions functions for Stokes / and Q for three representative cases, 
namely at 3871.2 A (bandhead), at 3869. 1 A (CN blend), and for 
the local continuum. 

The height dependence of the anisotropy at the bandhead 
wavelength clearly displays the role of the two competing ef- 
fects. Deep in the atmosphere (below km) the photons ther- 
malize and anisotropy vanishes even in presence of a large-scale 
temperature gradient. Around a height of 400 km the surface 
effect dominates and causes a strongly negative anisotropy. At 
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even higher layers, the surface effect diminishes and, as a conse- 
quence, the source function gradient effect manages to increases 
the anisotropy to low but positive values. The anisotropy remains 
low at the top because the temperature gradient and thus also the 
source function gradient are relatively small near the tempera- 
ture minimum. In ID modeling and for a given wavelength the 
anisotropy only varies with height. However, when observing 
closer to disk center we sample deeper parts of the atmosphere 
so that we get greater contributions from the layers with nega- 
tive anisotropy. Therefore, with increasing yU value, the value of 
the emergent polarization becomes smaller and even negative for 
high ^ (see Fig.|2]i. 

The behavior of the anisotropy for the CN blend at 3869.1 A 
is in principle similar to that in the bandhead, although the emit- 
ted radiation at this wavelength forms lower in the atmosphere, 
where the temperature gradient dT jdh is larger Also, the gra- 
dient in the optical depth scale dT/dr is even higher, because 
of the smaller opacity at this wavelength. It is this gradient that 
defines the limb-darkening law and, therefore, important for the 
anisotropy. Hence, the anisotropy at 3869.9 A is always greater 
than that at the bandhead and never becomes negative. 

The continuum radiation forms much lower in the atmo- 
sphere, where the temperature gradient is very high so that the 
surface effect becomes negligible. Therefore, the behavior of the 
continuum anisotropy is simpler than in the previous cases, and 
it monotonically increases with height. 

We considered a number of additional known models such as 
F ALA, FALC, FA LF and FALP of Fonten laetaL] ( |T993]l, FALX 
of |Avrett| ( |1995| l and colder models like model 2 of Anderson 
d 1989b and AYCOOL of |Ayres et aL] ( [T986l ) and |Solanki eTaL 
( 1994). Only the model AYCOOL which has a strong tempera- 
ture drop in the highest layers is able to produce relatively well 
the overall level of polarization, but it gives an incorrect polar- 
ization ratio in the bandhead compared to lines 0.5-1 A away 
from it. This model gives also completely wrong intensity pro- 
files. 

The main conclusion from the above discussion is that the 
anisotropy of the radiation field driven by the temperature gra- 
dient in the FALC model is too low to explain the observed po- 
larization. On the other hand, increasing the temperature gradi- 
ent, for example, by lowering the temperature in the chromo- 
sphere, fails to reproduce the observed center-to-limb variations 
in the intensity. Therefore, it appears that one-component mod- 
eling can reproduce either I/I^ or Q/I signals depending on the 
choice of the temperature gradient, but not both simultaneously. 

As was discussed in Sect. 13.41 the source function of the i- 
th line contains two scattering terms: the Rayleigh one, which 
corresponds to the radiative excitation by the photon absorption 
in the same /-th line and the Raman one, which corresponds to 
the radiative excitation by the photon absorption in another j-th 
line (which shares the upper state with the /-th line, see Fig.[Ti. 
Therefore the Stokes Q source function of the /-th line depends 
not only on the anisotropy of the radiation at the wavelength 
of the /-th line, but also on the anisotropy of the radiation at 
the wavelength of the j-th line. In some cases these lines can be 
up to about 50 A away from each other. It means that due to 
the Raman scattering the polarization signals at two relatively 
distant regions of the spectrum can be tightly connected with 
each other This significantly complicates the calculations as for 
computing the polarization even in the narrow spectral region 
one has to calculate the intensity and the anisotropy basically 
for the entire CN violet system. This system is blended by the 
large amount of the strong atomic lines, many of those have to 



be calculated in NLTE ( [Shapiro et al.|2010) l. For simplicity we 
calculated all the atomic lines in LTE but slightly adjusted their 
oscillator strengths. 

The effect from the Raman scattering is especially strong in 
the case of the optically thick lines which presence significantly 
modifies the anisotropy of the radiation field. For example with- 
out the Raman scattering term the polarization at the bandhead 
(where the optical depth is the highest) would be even lower 
than shown in Fig. [2 and negative for all /i-values. However 
this strong negative polarization is diluted with the polarization 
brought from the different regions of the spectrum by the Raman 
scattering. 



4.2. Two-component atmosphere modeling 

Next we try to fit the observations with a multi-component at- 
mosphere model. There have been many theoretical and obser- 
vational evidences for a temperature bifurcation in the upper 
photosphere and lower chromosphere. One comes from obser- 
vations of CO fundamental vibration-rotation transitions in the 
near infrared ( |Noyes & HaII|[T972| |Uitenbroek|[2000| and ref- 
erences therein). Also, to explain the difference between field 
strengths deduc ed from Hanle effect analys is of the Sr i 4607 



A and C2 lines, Trujillo Bueno et al. (2004) suggested that the 



observed polarization in these lines is formed in regions with 
different temperature and anisotropy (see also Asensio Ramos 



|& Trujillo Bueno|2005 1 [Trujillo Bueno & Shchukina|2007| and 



references therein). Further, Holzreuter et al.[ ( [2006| l showed that 
a temperature bifurcation is necessary for the explanation of the 
linearly polarized Ca 11 K profile. 

We tried to fit the polarization signal with the combination 
of the models listed in the Sect. 14.11 However the differences 
between the anisotropics in the considered models are very small 
(except AYCOOL model). The anisotropy have also a different 
sign in different parts of the solar atmosphere. We found that the 
combination of these atmosphere models also cannot reproduce 
the observed Q/I curve. 

We tried also to construct new two-component atmosphere 
models. We found that it is possible to construct a model which 
can fit reasonably well all five Q/I curves but at the same time 
fails to fit ///c curves. We found another model which can fit si- 
multaneously ///c and 2/^ curves for one particular yu-value, but 
fails to fit observations for other /i-values. So, we can conclude 
that although some compromises can be found, it is impossible 
to fit simultaneously the whole set of the observed data by com- 
bining two models. 

The failure to fit simultaneously I/I^ and Q/I signals with 
the considered combinations of atmosphere models for the broad 
range of //-values is not surprising. To reproduce the polarization 
signal the temperature gradient in one of the atmosphere compo- 
nents has to be larger than in FALC (see Sect. 4.1 1. However, 



such a large gradient will result in wrong intensity center-to- 
limb variations, so the temperature gradient in the second atmo- 
sphere component has to be lower than in FALC to compensate 
it, which in turn cause a very small or even negative polariza- 
tion signal. So the model which has to correct center-to-limb 
variations in the intensity spoils the polarization center-to-limb 
variations. Introducing atmosphere models with more than two 
components will not change the situation. 
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Fig. 4. The dependence of the anisotropy on the height calculated with the FALC model (solid lines) for several / parameter values 
(marked near the curves), and with the FAL A (dotted line) and FALF (dashed line) models for f - (this implies calculations, 
without anisotropy alteration, see Sect. 5.2 1. The anisotropy was computed at 3871.2 A in the (1,1) bandhead (left panel) and at 
3869.1 A in a CN blend (right panel). 



5. Effect of horizontal temperature variations 

5.1. Additional anisotropy due to 3D effects 

The main problem in the interpretation of our observations is a 
lack of radiation field anisotropy. There are many reasons for 
underestimating the anisotropy in the ID modeling. Basically, 
each heterogeneity in the solar atmosphere will increase the 
anisotropy. Using snapshots from a hydrodynamical simulation 
of the solar surface convection Shchukina & Trujillo Bueno 
( 2003| l found indeed very strong horizontal fluctuations of the 
anisotropy in the Sri line. It varies from negative values up to 
about 15% at the height of 400 km. 

One mechanism for increasing the anisotropy can be due to 
horizontal temperature variations. Let us consider, for example, 
a two-component atmosphere model. In general it is incorrect 
to separately calculate Stokes Q in each of the components and 
then linearly sum these individual contributions up to find the 
total emitted Stokes Q. The adjacent areas with diff'erent tem- 
peratures affect each other and change the anisotropy behavior 
Indeed, the areas of the atmosphere which are described by the 
warm component are surrounded by colder areas of the second 
cold component. Therefore, close to the boundary the radiation 
coming from aside originates in the cold atmosphere compo- 
nent, while radiation coming from the radial (vertical) directions 
forms in the warm component. As the anisotropy is basically 
a measure of the excess of the radial radiation over the radia- 
tion from the side, this will lead to an increase in the anisotropy 
as compared to the warm model. The same effect decreases 
the anisotropy in colder areas near the boundary. However the 
anisotropy alteration scale in some frequency v (which is a char- 
acteristic distance from the boundary to the regions where the 
considered effect can be neglected) is inversely proportional to 
the opacity at this frequency (as these scales correspond to the 
optical thickness one in both components). But the CN number 
density and, consequently, opacity strongly decrease in warm 



regions, so the increase in the anisotropy in such regions is 
significant over a greater geometrical scale than the decrease 
in cold regions. The considered effect will increase the mean 
anisotropy and therefore polarization as the anisotropy behavior 
at all depths has the same tendency: the total mean anisotropy 
is larger than the mean of the individual atmospheres. The il- 
lustration of the anisotropy behavior close to the boundary be- 
tween different atmospheric components is given in the available 
online Sect. |A] The efficiency of this effect in atomic lines and 
the continuum is much weaker, since their opacity dependence 
on temperature is not so steep as for molecules. Let us notice 
that the same effect changes also the mean intensity and, conse- 
quently, the intensity part of the source function and the emitted 
Stokes /. However, the relative impact on the mean intensity is 
much smaller than relative changes of the anisotropy (as for the 
mean intensity the sum of the intensities coming from the ra- 
dial and side directions is important, while for the anisotropy it 
is their difference). Moreover, the intensity source function con- 
tains also the thermal part (see Eq. (j7]l) while the polarization 
arises only in the scattering process. 

If the photon mean free path is much smaller than the char- 
acteristic size of alternating cold and warm areas, the considered 
boundary effect is not important, since the anisotropy in most 
of the area is equal to the asymptotic value. If on the contrary, 
the photon mean free path is big, the horizontal radiation comes 
from many alternating parts of the atmosphere and the effect is 
smoothed out. Therefore, the considered anisotropy variations 
are most efficient when the photon mean free path is comparable 
with the granulation scale. 

The mean photon free path in the layers with the strongest 
contribution to Stokes Q is about 300-500 km and, therefore, is 
of the same order of magnitude as the granulation scale. So, one 
can expect that the anisotropy alteration could be quite signif- 
icant in CN lines. However, in spite of the appropriate photon 
mean free path, this effect is not important for the continuum, as 
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it only weakly depends on temperature. Therefore, an increase 
in the anisotropy in the warm atmosphere is almost exactly com- 
pensated by its decrease in the cold component. 

Let us notice that the above consideration is valid only for 
optically thick lines. The CN lines under consideration remain 
optically thick even in the warmer atmosphere component. In the 
optically thin case the scattering polarization becomes propor- 
tional to the molecular number density and, therefore, the contri- 
bution from warm regions reduces. This situation was discussed 
by Trujillo Bueno et al. ( 2004} who found that the scattering po- 
larization in optically thin C2 lines mainly comes from colder 
upfl owing regions. 

The considered effect can be the reason of our failure to 
reproduce the observations with a two-component atmosphere 
model in Sect. 4.2 Such modeling underestimates the anisotropy 
and, therefore, scattering polarization. It is important to notice 
that this conclusion is independent on the real granulation model. 
It is not even essential for us which part of the solar atmosphere 
is warmer at the line formation height: granular upflow or in- 
tergranular downflow. The origin of the reversed granulation in 
the solar atmosphere and the place where it occurs is still not 
fully clear (see, for example Cheung et al.|2007 and references 
therein). 

It appears that a self-consistent 3D modeling can improve 
the situation with the anisotropy. But even current advanced 3D 
models may not be able to provide the correct anisotropy as 
they involve many free parameters and assumptions. Thus, their 
ability to reproduce the CN observations analyzed in this paper 
could be a test for their reliability. 

5.2. Anisotropy correction in 1D modeling 

To account for the underestimated anisotropy we artificially in- 
crease it by introducing a weighting function. Namely, in calcu- 
lations of the second moment we multiply the intensity with a 
function //(//) only within the Stokes Q component of the source 
function, 

+1 
-1 

+3(fi'^-l)Q(x',fi'))dfi', (26) 

which is chosen to overweight the radiation coming from the ra- 
dial (vertical) directions and give a lower weight to the radiation 
coming from aside. We used the following expression 



3(1 +V) 
3+/ ' 



(27) 



where / is a free parameter. Higher values of / imply 
stronger corrections, while /=0 corresponds to the uncorrected 
anisotropy. Such a choice of function is quite convenient as in 
the case of zero anisotropy our correction obeys photon conser- 
vation. In general the error caused by additional photons remains 
on the order of the anisotropy, i.e. a few percents and is thus 
negligible because we did not introduce any correction to the 
zero-order moment of intensity (see Eq.([T6|)). Note that when 
applying the correction factor riifi') also to Stokes / component 
it would be essential to normalize the scattering integral. We did 
not introduce any correction to the zero-order moment of inten- 
sity (see Eq. (16 1) and also kept the intensity fixed during the 
POLY iterations. With the purpose to avoid any additional dif- 
ferential effects we kept the parameter / fixed in the atmosphere 
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Fig. 6. -contours for solutions in the region of the (0,0) band- 
head. The ellipses correspond to discrete values of / indicated 
at the top. Three contours are drawn for each / corresponding to 
confidence levels 68.3%, 90%, and 99.73%. 



down to 150 km and then linearly decreased it to zero at -100 
km. Examples of the anisotropy behavior with depth for several 
/ values are given in Fig.|4]for the case of the FALC atmosphere 
model. One can see that already /=0.3 value can solve the prob- 
lem of negative anisotropy at the bandhead frequency, although 
as will be discussed in Sect. |6]for better fit we need a somewhat 
higher value. 



6. Spectral CN modeling with anisotropy correction 

In the previous section we showed that introducing the new pa- 
rameter / into Eq. pS] ) can solve the problem of the radiation 
field anisotropy deficiency. In this section we are employing 
again the standard FALC atmosphere model but with the / pa- 
rameter to regulate the anisotropy. Apart from this anisotropy pa- 
rameter, there are two other main free parameters: the fraction of 
thermal photons at the temperature minimum layer 6'™" and the 
magnetic field strength B. Variation of the first parameter affects 
all lines in more or less the same way, although some differential 
effects can be induced by different line formation depths and by 
the dependence of collisional rates on height. Effective Lande 
factors of the upper states for the CN violet system lines are 
inversely proportional to the angular momentum A^, and there- 
fore the magnetic sensitivity is decreasing with (see Eqs. (j5- 
|6]l). Thus, the Hanle eff'ect acts on CN lines differentially which 
greatly reduces the model dependence of the deduced field. For 
more details see the discussion in ,Shapiro et al.,(2007a) l. 

In our model the turbulent single-value magnetic field is as- 
sumed to be constant in the whole atmosphere. In Sect. 6.2 we 



will discuss the possibility to confirm or exclude some variabil- 
ity with the height. 

6.1. Region of fhe (0-0) bandliead 

Here we present our results for the spectral region close to the 
(0,0) bandhead (3880.5-3883.6 A). We found that for values 
/ > 0.4 the observation can be fitted with a good quality. 

In Fig. |5] we show one of the good solutions yielding a 
magnetic field strength of B = 82 G, a collisional coefficient 
<5jf|'" = 0.3 and an anisotropy parameter / - 0.6. One can see 
that both Q/I and I/I^ signals can be fitted well except for a few 
wavelengths affected by atomic blends. We recall here that back- 
ground atomic lines were calculated only in LTE. Comparing the 
observed and zero-field Q/I profiles one can see that the Hanle 
effect significantly depolarizes the spectrum at both fi values. 
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Fig. 5. Observations of Stokes ///c and Qjl in the region of the (0,0) bandhead and fits obtained with the FALC atmosphere model 
and increased anisotropy. The three curves presented in each panel correspond to observations (thick solid) and to the calculations 
with zero magnetic field (dotted) and with a magnetic field strength of 82 G (thin solid). Note that no differences between the 
calculated curves are apparent in intensity. 



Moreover its action leads to an apparent wavelength shift of four 
polarization peaks between 3880.6 A and 3881.8 A. The rea- 
son for this is that these peaks are formed by pairs of doublets 
with angular momentum about 50 and 10, respectively, and 
in each peak the doublet with smaller (higher Hanle sensitiv- 
ity) has a shorter wavelength and forms a blue part of the peak. 
Therefore, the blue part of each peak is more strongly affected 
by the Hanle effect than the red one leading to the apparent shift. 
Such differential line behavior in the magnetic field allows us 
to disentangle the polarization changes caused by the coUisional 
coefficient 5™" and magnetic field B and, hence, to determine 
these parameters. 

However, increasing the coefficient /, which is in charge of 
the additional anisotropy, will lead to an overall growth of po- 
larization. Such a growth can also be achieved by decreasing the 
collisional coefficient. Thus, it is difficult to distinguish these 
effects with observations at only two different limb angles as ad- 
justments of both parameters influence polarization in approxi- 



mately the same way (see, however, the discussion in Sect 6.2 



where we are dealing with observations at five different /i). The 
///c spectra are not so sensitive to the 5™'^ coefficient and, there- 
fore, cannot be used as an additional constraint. Finally, the ef- 
fect of the increased coefficient / can be compensated by in- 
creasing the collisional coefficient <5jf|'" and by increasing the 
magnetic field strength B. Therefore, we have executed each;^^- 
minimization with a fixed / parameter 

To estimate the quality of the fit we use the concept and 
define 



O-Q/I 



(28) 



where n is the number of data points. For constraining our free 
parameters we used only Stokes Q/I as intensity is not very sen- 
sitive to variations of / and 5™" and the line depths depend on 
the line blanketing parameters which can be slightly adjusted 



(see Sect. |2]i. To compute we used 120 data points chosen 
under the condition to avoid blending with atomic lines. 

The ;^'^-contours for different / parameter values are pre- 
sented in Fig. |6] For each / value we have drawn three con- 
tours which correspond to the confidence limits 68.3%, 90%, 
and 99.73%. Of course, the errors in our case are not necessarily 
distributed normally and there are definitely some systematic er- 
rors. The deduced magnetic field strengths corresponding to x^ 
minima are listed in Table [T] A minimum / parameter of 0.4 is 
required for simultaneous fits to I/I^ and Q/I. The minimum 
value corresponds to the magnetic field strength B = 82 ± 10 G. 

The main conclusion from Fig.|6]and Table[T]is that the mag- 
netic field strength appears to be model dependent. A larger ad- 
ditional anisotropy leads to a stronger field because the colli- 
sional rates enter into the denominator of Eq. (|6]l. It implies that 
their increase will lead to weaker Hanle sensitivity, so for ob- 
taining the same Hanle depolarization a higher field is needed. 
Moreover, the overall polarization growth can also be partially 
compensated by the Hanle depolarization (however, only within 
a limit which is constrained by the differential behavior of the 
Hanle effect). For a fixed / value both the magnetic field and 
collision rate are well constrained. Even though decreasing the 
collisional rate (and hence increasing scattering) can be partly 
compensated by the growth of Hanle depolarization, the mag- 
netic field strength is constrained due to the differential behav- 
ior of the Hanle effect. Knowing collisional rates more precisely 
would help to obtain a unique solution. 

6.2. Region of the (1-1) bandhead 

Here we present our results for the spectral region close to the 
(1,1) bandhead (3868.8-3871.5 A). A comparison of the ob- 
served and calculated spectra with a magnetic field strength 
B = 45 G, collisional coefficient = 0.16 and an anisotropy 
parameter / - 0.6 are shown in Fig. |7] This region contains 
many CN lines from both (1,1) and (0,0) bands and also sev- 
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Fig. 7. Observations (thick solid lines) of Stokes ///(, and QII in the region of the (1,1) bandhead and fits obtained with FALC 
atmosphere model and increased anisotropy. Two curves presented in each Qjl picture correspond to the calculations with zero 
magnetic field (dotted) and with magnetic field strength of 45 G (thin solid). 



eral strong atomic blends (among them two strong iron lines 
at about 3869.5 A) as well as one CH line at about 3870.1 A. 
Therefore, the overall quality of the fit is worse than that for the 
(0,0) bandhead region. The minimum values and correspond- 
ing magnetic field strengths are listed in Table[r| the -contours 
are shown in Fig.js] The minimum value corresponds to the 
magnetic field strength B = 45 + 15 G. 



The observations with higher yu values sample deeper layers 
of the solar atmosphere where the coUisional rates and, conse- 
quently, the denominator of Eq. (|6]l are larger. Therefore, these 
lines are less affected by the Hanle effect and less depolarized. 
This can be clearly seen in Fig. |7] for fi - Q.\ the Qjl curve 
is strongly depolarized in presence of a magnetic field strength 
of 45 G, while for ^ - 0.5 it is almost unaffected by the Hanle 
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Fig. 8. The same as Fig.lijfor the (1,1) bandhead. 



effect. The height dependence of the coUisional rates and the de- 
duced magnetic field are strongly coupled with each other For 
example a steeper collision rate dependence on temperature (and 
hence height) will lead to magnetic field strength, decreasing 
with height. Thus taken the current uncertainties in coUisional 
rates and radiation field anisotropy it is not possible to unam- 
biguously determine the depth dependence of the magnetic field 
strength. Improvements in collision theory, 3D modeling of the 
radiation field anisotropy, and observations at a greater number 
of limb positions would be required for better constraints of the 
model atmosphere and the free parameters. 

Interestingly, there is one more mechanism which constrains 
the deduced magnetic field strength for this region. The alter- 
ation of the coefficient affects Qjl curves at all ji values in 
approximately the same way. The magnetic field however mainly 
depolarizes lines observed at small /i. This differential behavior 
significantly helps to constrain both the magnetic field strength 
and the coUisional coefficient (assuming that we employ the cor- 
rect dependency of the coUisional rates on depth). In the (1,1) 
bandhead region the considered mechanism is even more impor- 
tant than the differential Hanle effect as this region contains a 
crowded mixture of lines with different J numbers, and the dif- 
ferential line behavior in the magnetic field is not clearly seen. 
This mechanism is illustrated in Fig. [9] where the dependence of 
the mean polarization on /i is plotted for several values of the 
coUisional coefficient 5™" and magnetic field strength B. The 
observed five points best coincide with the curve which corre- 
sponds to B = 45 G and 5™" = 0.3. 

From Figs.|6]and|8]and Table[T]one can see that the magnetic 
field strength deduced with the same / parameter for both spec- 
tral regions differ by a factor of two, with the one for the (0,0) 
bandhead being larger This can be a signature of spatial mag- 
netic field variations as the observations of these spectral regions 
were made on different days and sample accordingly different 
regions on the solar surface. We also interpreted the data from 
the atlas by Gandorfer ( 2005| l for the (0,0) bandhead region and 
found a weaker magnetic field of about 40 G compared to the 82 
G obtained from our new observations with / = 0.6. However, 
since the applied constraints are sensitive to the amount of avail- 
able data (e.g., limb angles) only simultaneous measurements at 
the same limb distances in both regions can clarify the situation. 



7. Conclusions 

We have employed a new computational scheme to interpret si- 
multaneously center-to-limb variations of the intensity and linear 
polarization in the CN violet system. An analysis of these vari- 
ations improves our understanding of the scattering polarization 




Fig. 9. The dependence of the mean polarization in the (1,1) 
bandhead region on the /i value for several coUisional coeffi- 
cients 5™" = 0.05 (soHd lines), 0.3 (dashed Hnes), 0.7 (dashed 
dotted lines), and magnetic field strengths: G, 20 G, 45 G, 60 G, 
and 90 G (from top to bottom) for each given value of 5™". The 
five stars correspond to the five observed mean values. 

Table 1. -minimum values and the least squares fits to the 
magnetic field strength B for several / parameter values. 



(0,0) band region 



(1,1) band region 



/ 


^ mm 


H,G 
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3.12 




5.03 




0.4 


1.12 


59± 10 


2.26 


30± 15 


0.6 


0.98 


82 ± 10 


2.06 


45± 15 


0.8 


1.07 


102 ± 15 


2.25 


57± 15 


1.0 


1.24 


119 ±30 


2.50 


70± 15 


1.2 


1.42 


133 ±30 


2.72 


82± 15 


1.4 


1.61 


146 ± 30 


2.91 


105 ± 15 



in the CN violet system. In our model we solve the statistical 
equilibrium equations and self-consistently account for multiple 
scattering in optically thick molecular lines. The use of the Hanle 
effect allows us to probe unresolved entangled magnetic fields. 

We found that the standard FALC model is able to correctly 
reproduce intensity center-to-limb variations. However it fails to 
fit the linear polarization, in particular close to the CN bandhead 
where the calculated polarization becomes even negative (per- 
pendicular to the solar limb). Moreover, we have shown that si- 
multaneous fits to the intensity and linear polarization center-to- 
limb variations are not possible in ID modeling, neither for one- 
nor for multi-component atmosphere models. Models which can 
provide the correct intensity center-to-limb variations are not 
able to produce enough anisotropy of the radiation field to fit 
the linear polarization. 
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We discuss several physical mechanisms which can increase 
the anisotropy. One of them is due to horizontal temperature fluc- 
tuations in the atmosphere. This is especiaUy valuable for molec- 
ular lines due to the strong dependency of the molecular number 
density on temperature. 

We account for the imderestimated anisotropy in ID mod- 
eling by introducing an empirical weighting function into the 
radiative transfer. With this method we are able to fit center-to- 
limb variations of the intensity and polarization simultaneously. 
The magnetic field can be determined via the differential Hanle 
efi"ect but its value depends on the appUed anisotropy correction. 

The uncertainty in the fitted parameters such as the mag- 
netic field is strongly connected to simplifications and shortcom- 
ings in the currently available theory of molecular collisions. 
Increasing the anisotropy in ID models requires higher depolar- 
ization, which will raise the low Umits of collisions and magnetic 
fields. The same efifect can occur in 3D models, if the anisotropy 
is overestimated, or vice versa. 

Scattering polarization provides additional and independent 
observational constraints for testing the vahdity of model atmo- 
spheres, not only in the case of ID models but in particular also 
for testing results of 3D MHD simulations. 
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Appendix A: Anisotropy variations 

An example of anisotropy variations at the CN bandhead fre- 
quency due to horizontal temperature variations is given in 
Fig. A.l for a very simple, purely illustrative model. We con- 
sidered two adjacent homogeneous and isothermal atmospheres. 
One with temperature 4300 K (marked "cold") and another 
with temperature 4900 K (marked "warm"). Such a choice 
gives a realistic temperature difference at line formation heights 
(Uitenbroek 2000 ). The densities in both atmospheres were cho- 
sen under the condition that they are in the hydrostatical equi- 
librium with the FALC model at the height of 500 km (where 
the bandhead radiation at// = 0.1 is formed). The anisotropy be- 
havior is drawn for four horizontal planes: at the surface (depth 
D = km), at the distance of one photon mean free path from 
the surface in the cold atmosphere {D - 200 km), at the distance 
of one photon mean free path from the surface in the warm atmo- 
sphere (D = 2000 km), and deep in the atmosphere (D - 10000 
km). We indicated in boldface the anisotropy behavior in the nar- 
rower region of 1000 km as this roughly corresponds to the gran- 
ulation scale. Far from the boundary the anisotropy asymptoti- 
cally reaches the limits individually calculated for each atmo- 
sphere (shown with dashed lines). At the surface, the asymptotic 
values are zero for both atmospheres (as the temperature gra- 
dient is zero), but the actual anisotropy is positive in the warm 
atmosphere and negative in the cold one. The behavior of the 
anisotropy in the warm atmosphere is opposite to that in the cold 
one, while the mean anisotropy is positive. At the depth of 200 
km both asymptotic limits become negative because of the sur- 
face effect. At the depth of 2000 km the asymptotic value for the 
cold atmosphere increases again to zero, as the layer above is 
optically thick and the surface effect can be neglected, while for 
the warm atmosphere the asymptotic value is still negative. And 
finally for the depth of 10000 km both asymptotic values reach 
again zero. 
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Fig. A.l. The dependence of the anisotropy on the distance from 
the boundary between cold and warm atmospheres. In each panel 
the cold atmosphere is shown to the left, the warm atmosphere 
is shown to the right. D indicates the depth from the top of the 
semi-infinite homogeneous atmospheres. 



